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Abstract 

We study the behavior of dynamical critical exponents of the two- 
dimensional Ising model with a line of defects. Simulations done at 
early time (first hundred Monte Carlo steps) reveal that the critical 
exponent 6 of Janssen et al depends on the strength of the exchange 
coupling constant(J') of the altered line. On the other hand, our sim- 
ulations permit to conclude that the dynamical critical exponent z is 
not sensitive to changes in J'. In addition, we investigate the possible 
invariance of the anomalous dimension (xo) of the magnetisation at 
the begining of the process. 
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A lot of results have been recently obtained concerning the critical dynam- 
ical behavior of statistical models at early time. This kind of investigation 
was motivated by analytical and numerical results contained in the papers 
of Janssen et al. |J and Huse Q pointing to a new dynamical critical expo- 
nent (9) describing the now called "critical initial slip" phenomena. That 
exponent describes the raising of the magnetisation at very short time 

M(t) ~ m</ (1) 

when a system initially in random states, with a small magnetisation, is 
quenched to the critical temperature and evolves with the dynamics of model 
A. Simulations done for the Ising and Potts (q=3) model indicate that the 
new exponent 9 does not depend on the kind of dynamics (Metropolis, heat- 
bath or Glauber). In addition, generalized Binder's cummulant || as well 
as other sample averages 0'0 which scale as L° have proved to be useful 
in determining the dynamical exponent z which relates time and spatial 
correlation length (r oc £ 2 ). The power of this approach was also verified in 
nonequilibrium models such as the majority vote model 0. In this letter we 
investigate the behaviour of the dynamical critical exponents [z and 9) when 
a marginal operator is present in the Hamiltonian. Although the presence 
of a marginal operator has well known consequences for the static critical 
behaviour |7j], very little is known for the dynamical one. To the best of 
our knowledge the only study on this subject was recently conducted by Li 
et al. H for the Ashkin- Teller model which interpolates between Ising and 
q = 4 Potts model, exhibiting a continuous varying exponent (u) for the 
correlation length (£). They have also detected a nonuniversal behaviour of 
the dynamical exponent z. In the present paper we study the two-dimensional 
Ising model with a line of defects (diferent coupling constant between spins 
along a given line) to verify the influence of the marginal operator on the 
dynamical exponents 9 and z. In fact, we found that along the altered line 
the exponent 9 depends on the ratio (J'/ J). On the contrary, we detected 
no changes related to the exponent z, a hypothesis that we used to obtain 
the dependence of the critical exponent rj on the ratio of couplings (J' /J) 
is. Finally, we have investigated the behavior of the anomalous dimension 
x of the magnetisation at the begining of the process. At least within the 
precision of our calculations, that value was kept constant no matter the 
value of the coupling constant (J'). In the sequence, we present the model 
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and the results of the simulations we performed. 

The two-dimensional isotropic Ising model undergoes a spontaneous sym- 
metry breaking || at a critical value of the temperature given by T = 
2J//cln2 . At the neighborhood of that temperature, the correlation length 
diverges as (T — T c ) _1 which means that v — 1. In addition, the magnetiza- 
tion vanishes like (T — T c ) 1 ^ which leads to the critical exponent (5 — 1/8. At 
long distances the correlation energy-energy behaves, at the critical temper- 
ature, as r 
the relation: 

(J'), along just one line of the lattice (d = 1), fulfils the necessary conditon 
x op = d = 1 to obtain nonuniversality the presence in the hamiltonian 
of a operator which scales as r~ d . This, in fact, was the result obtained by 
Bariev [[H]] and McCoy and Perk|ll]], later rederived by Peschel and Schotte 
in the context of quantum chains. After the work of Cardy ]T3 , point- 
ing the importance of the conformal invariance in two dimensional systems, 
the Ising model with a line of defects was revisited by Turban ||14| . In the 
sequence, many body techniques |15j were used to investigate its quantum 
analog [16|'[l7|' ||18|| but after all, this model continues deserving attention 



2Xe , where x e is the anomalous dimension of the energy satisfying 
x e + 1/v = 2. So, the introduction of a different coupling 



[!£}. It is interesting to notice that the "global" exponent v as well as the 
critical temperature does not change when the value of the coupling (J') is 
altered. The same is not true for the exponent [jf = 2(3*) of the correla- 
tion function along the line of defects (at T = T c ). This "local " exponent 
depends continously on the value of the coupling (J ') as 



r]* = {-tan^Kf 



(2) 



where 



K = tanh(J!/A;r c )/tanh(J/A;T c ) 

and Ji is the dual of the modified coupling J' 

exp(-2Ji/A;T c ) = ttmh(j'/kT). 

Our investigation began by following the evolution of the magnetisation 
when samples are sharply prepared with mo but small. We measured 
the average over samples (10 5 ) of the magnetisation for two lines (one pure 
, another with defects). In addition, we also calculated the magnetisation 
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over the entire sample. Simulations were done for lattices with 36 < L < 72 
and the update followed from heat-bath dynamics. Because our interest 
is in the developing of the magnetisation of the lines, we prepared all the 
samples with the same initial magnetisation for each line. In Figures l.a and 
l.b we present the raising of the magnetisation for that three cases: pure 
line, line with defects and entire sample, when J' /J = 2 and 0.5. We call 
the attention of the reader for the similarity between the curves associated 
to the pure line and to entire sample in both cases. Even working with 
small lattices we note that the effect of the line of defects restricts itself to 
its neighbourhood. The conclusion is that the exponent 9 depends on the 
coupling constant of the defect line (see Figure 2) as occurs with the static 
critical index 77*. Dynamical universality (same results with heat-bath and 
Metropolis updating) was observed in the simulations as found in the pure 
Ising model PH , although bigger deviations could be seen when J' / J ^ 1. 

Next step is to check the plausible hypothesis that the dynamical critical 
exponent z is not changed by the line of defects. As pointed by Li [plU , 
the short-time behaviour of the magnetisation, of samples that at t — are 
totally magnetized (m = 1), is given by 

m ~ (3) 

where z is the dynamical critical exponent, 77 and u the known static indices 
of the model. If our hypothesis is correct we can obtain the exponent 77* which 
characterizes the polynomial decay of the correlation function (at T = T c ) 
along the line of defects by comparing the slope of two straight lines in the 
log-log plot of magnetisation versus time (one of them belonging to the line 
of defects, the other corresponding to a normal line far enough of the altered 
one). In Figure 3 is shown the behaviour of the magnetisation of those 
lines, for several values of the coupling constant J', whereas in figure 4 are 
presented estimates for the exponent 77*, obtained by the relation 

tana' 

V = V7 4 

tan a 

where 77 = 1/4 is the critical exponent for the pure Ising model. To obtain 
that equation we have used that z has ever the same value, no matter the 
value of J' is. In addition, tana'and tana are the slopes of the curves for 
altered and pure lines. The results so obtained are in complete agreement 
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with the exact ones (continuous line) obtained by Bariev JTU| and McCoy 



and Perk [11]. To obtain the error bars (smaller than the size of points) 
we repeated five times each simulation. We observe that our estimates are 
independent of the specific value of the dynamical critical exponent z. 

Another way to get the same conclusion is to calculate the evolution of 
the second moment of the magnetisation of samples which at t — satisfy 
the conditions : m = and £ = 0. Scaling arguments |§ asserts that m 2 
should behave as 

< m 2 >~ t( d -"A0/* ( 5 ) 

where d is the dimension of the object for which we are calculating that 
quantity (d — 1 when we are interested in lines). By saving the results for 
several values of J'/ J we have obtained the Figure 5 which shows polyno- 
mial behavior with distinct exponents. Finally, we can estimate rj* from the 
equation 

. tan</>' . , 

V * = l-(1- V )-^- 6 
tan0 

where tan<// (tan0) are the slopes of the straight lines in the figure, corre- 
sponding to the line with (without) defects. Results for 77* obtained by this 
procedure are plotted in Figure 6. Once again we see that the hypothesis of 
constant z works well, at least when J' /J < 1. The disagreement observed 
for higher values of the ratio J'/ J we atribute to the influence of the strong 
coupling over the fluctuations of the magnetisation of the pure line. 

Finally, we check the behavior of the anomalous dimension x of the 
magnetisation far from equilibrium. As it is now known, the critical initial 
slip is a consequence of the difference between the anomalous dimension 
Xq and the usual scaling dimension r]/2v of the magnetisation [[2!J. More 
specifically, we have 



x = r]/2u + z6 (7) 

what permits us to obtain xq ,for any value of the ratio (J'/ J), since we 
know both: the exponent rj (Figure 4) and the new dynamical exponent 9 
(Figure 2). As shown in Figure 7, the result is very approximately equal to 
0.53 (corresponding to the pure case) for a large interval of (J'/ J). We have 



5 



used the value 2.172 ±0.006 ,obtained by Grassberger[23], for the dynamical 
critical exponent z. 

In summary, in order to understand the effect of a marginal operator over 
the dynamical critical exponents we have investigated the short-time criti- 
cal dynamics of the Ising model with a line of defects. We have confirmed 
nonuniversality for the "local" exponent 9 but not for the global dynami- 
cal exponent z. The universality of z permitted us to determine the static 
nonuniversal index 77*, of the correlation function along the line of defects, 
using just the hundred first steps of the simulation (free of critical slowing 
down phenomena). We have also detected an apparent universal behavior for 
the anomalous dimension xq of the magnetisation even at the line of defects. 
Although the numerical results can be a little bit different, they are qual- 
itatively the same when we use heat-bath or Metropolis updating. At the 



present, we are studying global and local persistence|]24j| ' ' ||26|| phenomena 
in this model. 
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Figure Captions 

Figure 1 - Log-log plot of the magnetisation versus time for (a) a pure line, 
(b) the line of defects, and (c) the complete square lattice. Samples (100000) 
were originally at high temperature but they had a small magnetisation per 
site m = 0. 0278. The size of the lattice in this case is L = 72 and the 
updating was done with heat-bath dynamics. Figure l.a shows results for 
the case (J' / J) = 2 whereas figure l.b presents results for (J'/ J) = 0.5. 

Figure 2 - Dependence of the new critical exponent 9 on the ratio of 
couplings (J'/J). The continuous curve is just a guide for the eyes. 

Figure 3 - Log-log plot of temporal decay of the magnetisation at the 
line of defects for several values of the ratio of couplings (J'/J). The slope 
of those curves depends on the ratio (J'/J) and give estimates for 77*. (Eq. 
(fD) Heat-bath dynamics was used to update samples (50000) which were 
completely magnetized at t = 0. The size of lattice was L = 36. 

Figure 4 - Estimates for the nonuniversal critical exponent if (charac- 
terizing the decay of the correlation function at the defect line) obtained by 
the scaling relation m ~ r"/ 2 " z , exhibited in the Fi gure 3 and Eq. (^). The 
continuous line is the exact result for r/*, given by Eq.@. 

Figure 5 - Log-log plot of the squared magnetisation at the defect line for 
several values of the ratio of couplings (J'/J). Scaling relations predict that 
the slope of the lines should be given by (d — 77* jv)jz. 

Figure 6 - Estimates for the nonuniversal critical exponent 77* obtained 
with Eq.(||). When J'/J is greater than 1 results are poor. 

Figure 7 - Plot of the anomalous dimension of the magnetisation xq = 
z6 + r]*/2u versus the ratio of couplings (J'/ J). 
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